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This annual Technical Report presents research conducted during 
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This study Is being performed within the Composite Materials 
Research Group at the University of Wyoming. The Principal Investigator 
Is Dr. Donald F. Adams, Professor of Mechanical Engineering. Mr. J. M. 
Mahlshl, Ph.D. student in Mechanical Engineering, performed the work 
contained herein, as part of his research as'^is. antshlp duties within 
the Composite Materials Research Group. 
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SECTION 1 


INTRODUCTION 

The present report Includes work performed during the third year 
of a NASA-Lawls grant to study the energy absorption mechanisms during 
crack propagation In metal matrix composites. 

The report of the first-year work [1] contains a literature review 
Cv'verlng the general area of mlcromechanlcs analyses of unidirectional 
composites, as related to the present study. During the first year, 
an existing elastoplastlc, generalized plane strain, finite element 
mlcromechanlcs analysis [2-4] was modified to Include crack propagation, 
follovrlng the general procedure developed earlier by Adams [5-8]. Both 
longitudinal and transverse cross section models were used to study the 
influence of a broken fiber on inelastic stress distributions. 

During the second-year study reported in Reference [9], the generalized 
plane strain crack propagation procedure was refined, and detailed results 
obtained. The analysis was also reformulated for an axlsymmetrlc model. 

This alternate two-dimensional analysis retained all of the general 
feature j of the generalized plane s:raln version. Its Intended primary 
application was to permit the study of a simple model composite consisting 
of a single broken fiber in a circular cylindrical sheath of matrix material. 
The boundary condition selected during a crack propagation Increment was 
constant applied stress. While a realistic condition. It resulted in a 
crack propagating catastroplcally once It Initiated. 

During the present third-year study, this boundary condition has been 
altered so that a constant boundary displacement can be maintained during 
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an Increment of propagation. Thug, tha applied atreaa reduces as the 
crack provagatas, resulting in crack arraatmant. This boundary condition 
Is representative of what can be readily aimulated In a laboratory 
environment, using a testing machine in its displacement-control mode. 

Having the ability to simulate the propagation of a stable crock, 
the present report focuses on the energy dissipation •’nd fracture response 
of a single broken boron fiber In a sheath of aluminum matrix material, as a 
composite model which could be readily correlated with experimental dnt.i. 


SECTION 2 


AXISYMMETRIC FINITE ELEMENT MODEL 

One purpose of the axlsynanetrlc model as deecrlbed in Reference [9] 
was to provide supporting correlations with the two-'dlmensional, generalized 
plane strain longitudinal and traverse cross section models. However, an 
even more Important purpose 's to provide analytical results for comparison 
with proposed experiments under similar conditions, i.e., a single broken 
boron fiber surrounded by a uniform annular sheath of aluminum matrix, as 
shown In Figure 1. Tlie formulation of the required ax^ symmetric element 
and tho corresponding computer program ore presented in Reference [9]. The 
finite element and the computer program were fully developed and some pre- 
liminary results were presented. However, it was noted that the initiation 
of a crack in the single fiber axlsymmetric model led almost Immediately to 
a catastroplc failure of the composite. Thus, the results were not useful 
from an experimental correlation point of view, which requires the deter- 
mination of measurable quantities such as crack opening displacements and 
surface strains during the process of crack extension. For experimental 
verification purposes, a constant boundary displacement crack propagation 
scheme has been Implemented In the program during the present study, which 
ensures stable crack propagation. 

It was also felt that the number of elements used in the previous 
study [9J was inadequate to represent the high stress gradients which occur 
near broken fiber ends. A refined model was evolved, shown in Figure 2, 
having 8A0 elemenls, nearly 400 of which are concentrated near the broken 
fiber end. An automated triangular element mesh generation routine [10] 
has been Implemented for more accurate and efficient grid data generation. 


Plane of 
Fiber Break 
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Figure 1. Axlsyrometr tc Analysis Model of a Single Broken Fiber In 
Sheath of Matrix Materia! 
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Figure 2. Finite Element Grid Utilized 
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A (ipcclal double node concept hae been used at the Junction of the 
broken fiber end and the surrounding natrlx» as shown In Figure 3» to more 
occurately represent actual physical conditions. A routine hus been 
implemented In the main program which automatlcaily regenerates the mesh 
data accounting for the double nodes. 

2.1 Constant Dlspl«*cement Loading Scheme 

The cataatroplc failures observed in the previous model f9] were 
attributed to the foot that the crack initiated by the first failure of 
an element led to failure of additional elements during the process of 
transferring the reaction loads of the failed element to the adjacent 
elements. Further failure of a very large number of elements was observed 
during the process of reaction load transfer of the failed elements, 
growing the crack r*pldly to the edge of the model at the same applied 
load which Initiated the fli'at failure. 

A simple scheme was developed during the present study which prevents 
such rarid growth of the crack. In this new scheme the loaded boundary of 
the specimen is constrained from further movemenl during the ensuing crack 
propagation after initiation. The boundary Is held at constant displacement 
until a state of equilibrium Is attained, l.e., a stress state at which the 
reaction loads of failed elements applied to adjacent elements do not lead 
to any further failure of elements. The applied loads are automatically 
continuously reduced during crack propagation to maintain eq.'lllbrlum. 

The procedure adopted in the computer program for the constant 
'.Isplacement condition Is as follows. After assembling the stiffnesses of 
each of the individual elements into a banded global stiffhess matrix of 
the form [2] 


|i^£^i?i;«i9^«?^fiPisifSiSi9 
^^i?k9iSiiS<Ci^lSllSllSlli9l9llSll9 
5i<i.'<i?^i?^lfllSieilfl8!llBI0 
i‘^k.'k^^t?£'t:«lSil9l9l|CilCllClllSl 
' o>c> p.'C' >:^>:^i^i^icii^isiis<i 

- ' ^ "vW*» I 
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Figure 3. Double Node SyHtom Uncd at the Broken Fiber End 
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where (F) ie the load vector, (u) It. the displacement vector, and [K] 1 b 
the globol BtlffncsB matrix. 

The maintenance of a condition of uniform axial displacement of both 
fiber and matrix at the loading boundary necessitates the use of the Branca 
boundary condition technique [11], This technique as adapted to the 
mlcromechanlcs analysis has been explained In detail In References [I, 2, 9]. 
In this technique, the equations representing the axial displacements of all 
loaded boundary nodes are added to the equation corresponding to the axial 
displacement at the corner of the region of Interest (l.e., to che axial 
displacement equation representing the extreme right loading edge boundary 
node, which Is numbered last for convenience In Identifying It). 

The constant boundary displacement condition Is easily applied by 

setting the value of the stiffness in the last equation, which corresponds 

to the displacement of the entire loading boundary, to a very large number 

30 

after the Initiation of the crack. A value of 10 psl has been used in 
the present analysis. 

The value of the load reduction which results when the displacement L;’ 
constrained Is obtained by multiplying the displacement vector with the 
stiffness coefficients corresponding to the last boundary node axial 
displacement, l.e., 

{? } - [kKu) 
z 


( 2 ) 
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2.2 Boundary C'^nditlons 

The constant lateral dlsplaeamant boundary conditions applied In the 
original version of the axlsymmetric model (9] have been removed as they 
do not adequately represent physical conditions. The outer surface of the 
matrix la now free to deform in any direction. 

During the present analysts, It was found that the common node point 
for both fiber and matrix at the outer radius of the broken fiber, on the plane 
of the fiber break, led to very high stress concentrations in the fiber 
material. The reason for the presence of these larg« stress concentrations 
was found to be due to the modeled continuity of the vsupposedly completely 
broken fiber at the common node. In order to represent the actual conditions 
of total discontlnui ty of the fiber at the break while retaining the 
continuity of the matrix material at the same point, a double, node approach 
has been incorporated. Two separate node numbers are assigned at the same 
point (Figure 2), one being considered to be aaso<^lated 'th the fiber 
material and the other with the matrix material. The boundary conditions 
are applied such that the node point In the fiber is allowed to have both 
axial and radial displacements (i.e., it Ls free to move In noth tl^e axial 
and radial directions) while the node point In the matrix is constrained 
by the symmetry conditions to have only a radial displacement. Actually, both 
of these nodes (Nodes 8 and 9 in Figure 2) should also be constrained to 
have the same radial displacement . However, this would lead to the task of 
implementing a constraii ed equation canabllity in the program. By experimenting 
with the program it was found that the affect of the absence of this constrained 
displacement condition on the resul ts was not significant and could be neglected. 
In practice, the element In the matrix at the br' ,.»n fiber end typically falls 
early In the loading process, making this node Ineffective as the element itself 
will not be effect::ve any lonper. 
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2.3 Failure Criteria 

The failure of Interest In the present study Is that of the matrix 
material, which Is exhibited as a matrix crack Initiation at the tip of 
the pre-existing crack or fiber break. This crack then propagates with 
Increasing applied load, corresponding to successive failures of finite 
elements in the vicinity of the crack t 

The matrix Is typically an Isotropic material (although the analysis 
Is not restricted to this assumption). There are a considerable number 
of failure criteria available, many of which have been used In connection 
with ccmposlte material analyses. 

In the prior version of the present axlsymmetrlc mode] analysis tl» 9], 
an octahedral shear stress failure criterion was used. In the present 
study, it was found that this failure criterion did not always lead to 
realistic results. Thus, the analysis and related computei program was 
modified to permit the selection of any one of six different failure 
criteria, at the user ' s option. These six criteria are presented In Table 1. 
Example results and additional discussion will be presented later, In 
Section 5 of this report. 

The use of the axlsymmetrlc model, with all of these modlilcations 
Incorporated, Is demonstrated by a series of numerical examples In the 
following sections. In the next section, SjcLion 3, the properties of the 
boron fiber and aluminum matrix utilized are presented. These are the same 
as those used in the two prior reports [1,9]. Thus, direct comparisons can 
be made for all results presented here. 

The axlsymmetrlc finite element model is first applied to a solid aluminum 
rod containing a central penny-shaped crack, in Section A.l. The results are 
compared with those obtained using other (closed form) methods, thus esta- 
blishing the performance of the finite element model. Then, In Section 4.2, 
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Failure Theory Options Available in the Present Analysis 



X, V, S are ultimate axial, transverse and shear strengths, respectively 
Ultimate octahderal shear stress (/I'/3 X ) 
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detailed resulte are presented for various boron /aluminum composite 
models. Finally, in Section 4.3, energy absorption associated with crack 
propagation Is discussed. 


SECTION 3 


MATERIAL PROPERTIES 

In modeling the boron/aluminum composite, the boron fibers have 
been treated as brittle, linearly elastic materials with isotropic 
strength and stiffness properties. The aluminum matrix has also been 
considered to be Isotropic, but Is modeled as an elastoplastlc material. 

To accomplish this, the actual stress-strain curve of the aluminum alloy 
selected Is Input to the analysis by curve fitting via a Rlchard-Blacklock 
two-parameter equation [21], as discussed in Appendix A-5 of Reference [1]. 
Thus, at any load level the tangent modulus for any given element can be 
computed. This makes possible an accurate representation of the plastic 
deformation of the matrix. 

Although the nonlinear material properties of any matrix material, 
e.g., another aluminum alloy, can readily be Incorporated in the analysis, 
a 6061-T6 aluminum alloy at 75“F was used In obtaining the present results. 
The material properties shown In Table 1 were obtained from Reference [22]; 
the full range stress-strain curve for determining the curve-fit parameters 
used Is shown In Figure 4. 


Table 2 

6061-T6 Aluminum Alloy Matrix Material Pvipertles [22] 

Young's Modulus E * 10.0 x 10^ psi 

Poisson's Ratio v * 0.33 

Tensile Yield Strength F^Y ■ 36000 psi 

Tensile Ultimate Strength F^^ * 45000 psi 

Coefficient of Thermal Expansion a * 13.0 x 10"^/®F 










Stress, (7 (ksi) 



The boron fiber properties Indicated In Table 3 were obtained from 
Reference [23]. 


Table 3 


Boron Fiber Material Properties [23] 

Young's Modulus E ■ 60.5 x 10^ psl 

Poisson's Ratio v *0.13 

Tensile Ultimate Strength Ftu . pty - 500,000 psl 

Ultimate Strain e^u « . g,264 x 10“3 in. /in. 

® .ft 

Coefficient of Thermal Expansion a ■ 9.0 x 1C "/"F 



Figure 4. 


Typical Full Range Stress-Strain Curve for 6061-T6 Aluminum 
Alloy at Room Temperature [22]. 


SECTION 4 


NUMERICAL RESULTS 

4.1 Evaluation of the Axlsymmetrlc Finite Element Model for Crack Problems 

In General 

According to linear elastic fracture mechanics (LEFM) theories, a crack 
will begin to propagate when the amplitude of the stress field In the 
immediate vicinity of the crack tip (the stress Intensity factor) reaches 
a critical value. Thus, the prediction of stress Intensity factors for 
different geometric conditions assumes major Importance In LEFM theories. 

Even though the present finite element model and related computer 
progr?im are not limited by linear elasticity assumptions. It was decided 
to first evaluate the performance of the model under this most severe case 
of singular stresses near the crack tip. 

4.1.1 Application of Conventional Elements 

A considerable! amount of work has been done in the area of finite 
element appllcatlDns to LEFM problems, and a number of two-dimensional and 
three-dimensional elements with embedded singularities are now available. 
However, conventional elements, even though they do not adequately represent 
the singular state of stress, can be used for rather accurate estimations of 
stress intensity factors by proper interpretation of the results. 

There are two different techniques for calculating stress intensity 
factors from conventional element solutions. The first technique requires a 
very fine grid near the crack tip to obtain a very accurate description of 
the displacement field. These displacements are then substituted into the 
classical continuum solution to obtain the stress Intensity factors. In 
the second technique, the strain energy of the system with a crack is 
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calculated for two slightly different crack lengths. The strain energy 
release rate during the crack extension gives the stress Intensity factors 
directly. 

4.1.2 Axlsywmetrlc Crack Problems 

The two types of axl symmetric cracks In an axl symmetric body that can 
be analyzed using the axlsymmetrlc finite element method are a central 
penny-shaped crack and a circumferential crack. In the present verification 
example, only a central penny-shaped crack In an aluminum rod subjected to 
axial tension will be utilized. 

The basic equations of linear fracture mechanics are [11,24] 

0^ - V cos((j)/2)[l - sln((|)/2)sin(3(j./2) J 

(2p)® 

o - ^ cos((J'/2)fl + sln(<(>/2)8ln(3(<./2)] 

^ (2p)* 

T * Y sln(i(i/2)c08(())/2)cos(3(|)/2) (3) 

(2p) 

K (2p)^ 

— [(5 - 8v)co8((Jl/2)-cos(3(^/2) ] + ^(l-v)cos (j) + . . . 

K.(2pl^ 

“z * — 8G — f (7-8v)sln((fi/2)-sln(3(f/2) I - sin .f. + . . . 

where 

Kj ■ stress Intensity factor 
p , (<i “ polar coord ..nates , as shown In Figure 5 
G “ shear modulus, v * Polrson's ratio 

a “ coefficient of second term in theasymptotlc expansion [25] 

0 , o , T * radial, axial and shear stresses 

rr zz rz 

u^, * radial and axial displacement components 
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The BCresH Inteneliy factor Kj Is calculated from the axpresslons for 

u and u by substituting corresponding values of u and u obtained from 
r 5! r z 

the finite element )Olutlon. 

For the above case of a central penny-shaped crack In an Isotropic 
thick cylinder subjected to Mode I type loadlngt there Is also a closed 
form solution available for stress Intensity factor, as follows [26-28]: 


K, 


a/va 

l-(a/b)' 


r— 

2 

1 + ^ - 

5a^ ' 

3~ 

4- 0.268 ~ 

TT 

2b 

8b^ J 

b^ - 


(4) 


where a, b and c are geometric parameters, as shown In Figure 5. 


A . 1 . 3 Stress Intensity Factor s 

Penny-shaped cracks of three different radii were analyzed. The 
overall finite element grid geometry was maintained; the radius of the 
crack was varied by changing only the boundary conditions. No effort was 
made to concentrate more elements near the crack tip as the purpose of 
this linear elastic fracture mechanics example was not so much to 
accurately determine the stress intensity factor, but rather to evaluate 
the finite element model to be subsequently used for studying composite 
behavior. Even better estimates of stress intensity factors would have 
been possible If that had been the primary purpose. 

The estimated stress Intensity factors calculated using the present 
finite element model, and the corresponalng values obtained using the 
closed form solution of Eq. (4), are presented In Table 4. 
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Table 4 

StroBB Intensity Factors 


ra/tb 

(see Fig, 5) 

K,/» 1 

present finite 
element analysis 

closed form solution 
(Eq. 4) _ 

0.25 

0.531 

0.819 

0.46 

0.809 

1.112 

0.84 

1.647 

1.698 


Since the stress Intensity factors are calculated using the displace- 
ment ccmponants obtained from the finite element solution, they are very 
sensitive to the accuracy of the displacement field. Since the constant 
st’.aln triangular element only allows linear variations of displacements 
within each element, a better representation of the displacement field 
can be achieved by having a larger number of smaller elements near the 
crack tip. 

s • 1 • 4 Inelastic Crack Propagation In an Aluminum Bar 

Having verified the accuracy of the axlsymmetrlc finite element analysis 
and related computer program In predicting the high localized stress concen- 
trations at the crack tip of a homogeneous, linearly elastic material, as 
presented in Section 4.1.3, the capability of the analysis to model inelastic 
material response and crack propagation was exercised. All of the modifi- 
cations described In Section 2 were incorporated. 

The example of the penny-shaped crack in an aluminum bar of circular 
cross section, used In Sections 4.1.1 through 4.1.3, was extended by 
increasing the applied axial stress to cause Inelastic material response 
and then crack propagation to failure. The stress-strain response predicted 
Is presented in Figure 6. It should be noted that the stress plotted Is the 
axial stress applied at a large distance from the crack site, l.e., the total 


APPLIE 



AVERAGE AXIAL STRAIN. SjCIO'^) 


Figure 6. Stress-Strnln Response of an Aluminum Circular Rod with 
Central Penny-Shaped Crack 
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applied load divided by the groea croea-aectlonal area. The first abrupt 
drop In applied axial streaH corresponds to the Initiation of crack growth* 
which becomes arrested as the applied stress drops during the constant 
boundary displacement crack grrwth Increment. The applied stress is then 
subsequently Increased as additional axial loading increments are applied* 
leading to additional Increments of crack growth* until the crack propa- 
gates completely outward across the entire cros^i section. 

4.2 Crack Propagation In a Bo rorV Aluminum Model Composite 

The constituent material properties presented In Section 3 were utilized. 
The single broken fiber model Is shovm In Figure 1, and the finite element 
grid utilized Is presented in Figures 2 and 3. 

The variable studied was the thickness of the aluminum sheath surrounding 
the broken boron fiber. This has been defined* as in the previous report [9], 
as the ratio of the radius of the fiber to the radius of the matrix sheath* 
l.e., rf/r,n. Three arbitrary radius ratios were analyzed, viz, rf/tm ■ 0.25, 
0.46, and 0.84, which correspond to fiber volume contents of 6.25, 21.2, 
and 70.6 percent, respectively. 

As In the prior work, no attempt was made to model thermal stress 
effects due tc fabrication processes, although the analysis has this 
capability. The maximum normal stress failure criterion was used as governing 
the crack propagation process (see Table 1). 

Results are presentea in Figure 7. As can be seen, there is a considerable 
amount of stable crack propagation, as indicated by the many abrupt drops in 
the applied axial stress. This increased with increasing thickness of the 
aluminum matrix sheath, as would be expected. Since the area under the 
stress-streia curve is proportional to the energy absorption capacity of the 
composite, this stable crack growth is obviously beneficial. 
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Figure 7 
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StrcHH-Strain Response of a Boron/Al umlnum Composite with 
Broken Fiber Under an Axial Applied Stress 
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Plct\')»'lal representations of the crack propagation patterns tor 

each of the -,/r_^ ratios of Figure 7 are shovm In Figure 8. The failed 
t m 

elements are shown shaded. The two values of applied axial stress (In 
units of psl) given under each diagram Indicate the stress level (from 
Figure 7) at which that Increment of crack propagation initiated, and 
arrested. As can be seen, the cracks tended to grow radially outward, with 
a slight tendency to also grow along the fiber-matrix Interface. 

The predicted ultimate axial st'ength for each of the rf/r„, ratios 
can also be obtained from Figure 7. These three values are plotted in 
Figure 9. The value for “ 0 represents the ultimate strength of 

the aluminum matrix (Figure 4 and Table 2), viz, 45 ksl. The value for 
“ 1.0 of zero corresponds to the trivial case of a broken fiber 
and no aluminum matrix. 

While a linear relation between ultimate strength and rf/r^, Is not 

necessary fot Jils nonlinear material response example, the data of 

Figure 9 do Indicate a relatively linear relation. 

An Important, longer term objective of the present study Is to provide 

numerical results which can be used to correlate with experimental work as 

It becomes available. Crack opening displacement would be one such 

experimental measurement (using X-ray or some similar technique for ..he 

present case of the optically opaque aluminum matrix). Figure 10 Is a 

plot of crack opening displacement (defined here as the distance between 

the two ends of the fiber at the break), for all three r./r ratios. As can 

t m 

be seen, uhe predicted displacement values are large enough (in thousandths 
of an inch) to be readily measured experimentally. 
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It would b« simpler, of course, to measure the surface displacement 

(strain) of the broken-fiber model at the plane of the fiber break. This 

meusurement could be made directly, by using a strain gage or extensometer. 

Predicted axial strains, e , on the surface of the aluminum matrix at the 

z 

plane of the fiber break, are plotted In Figure 11 for all three 
ratios. The corresponding circumferential strains, e^, are presented In 
Figure 12. The straight lines corresponding to the b train response of 
a solid aluminum rod are also given, for comparison purposes. 

As can be seen In Figures 11 and 12, the measurement of surface 
strains Is not a very sensitive method of determining the extent of crack 
propagation. 

A practlc^il consideration In attempting correlations with experimental 
data is how long the numerical model must be to avoid influences of end 
effects. It is usually deslreable tvO minimize the length of the model, 
in order to reduce the number of finite elements required, and thus to 
reduce computing time. The length of the model required can be determined 
by selecting a length, and then observing how rapidly the stresses return 
to their undisturbed values away from the discontinuity (the fiber break 
In the present case). The models used In the present examples were all 
4.1 composite model diameters long. 

In the case of a single broken fiber In a matrix sheath, the normal 
stress In the axial direction on a transverse cross section of the fiber 
remote from the discontinuity will be uniform, as will the norimil stress 
In the matrix sheath. The shear stress will be zero, including the shear 
stress at the fiber-matrix Interface. 

% 

Figure 13 Is a plot of the axial normal stress distribution across 
both the fiber and matrix, for the case rf/r^, ■ 0.46, on cross sections 
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Figure 11. Surface Axial Strain at Plane of Fiber Break Versus Applied 
Axial Stress 
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Figure 12. Surface Circumferential Strain at Plane of Fiber Break Versus i 
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at Increasing distances from tha fibar braak. Tha applied axial strass Is 
4 ksl, well within tha elastic range of material response. Right at the 
plane of the break, the axial normal stress In the fiber Is saro, by 
definition of n free surface. However, the stress concentration In the 
matrix adjacent to the fiber break Is severe, an axial normal stress of 
28 ksl being Indicated In Figure 13. However, at a distance of about 
only one and one-half fiber diameters, the stress distribution In both fiber 
and matrix Is almost uniform again. 

Figure 14 Is a plot of the shear stress along the fiber-matrix 
Interface, plotted as a function of distance In fiber diameters from the 
fiber break. It Is this shear stress which transfers the applied loading 
across the fiber break, via the matrix sheath. As can be seen, the shear 
stress builds up very rapidly near the break, but also decays rapidly 
with Increasing distance from the break. At three fiber diameters, the 
shear stress is negligibly small. 

The critical fiber length for complete load transfer, as predicted 
by the simple mechanics of materials models. Is 


9 , 


c 



1 1 l 

2 ^m 


(5) 


where Is the fiber axial strength and Is the matrix shear strength; 
as discussed by Chamls [29]. For the present boron/ aluminum composite, 

0 ^ ■ 500 ksl, ® 26 ksl, d^ ■ 0.0056 in. Thus, using Eq. (5), “ 9.6, 

corresponding to a distance of 4.8 fiber diameters from the fiber break In 
Figures 13 and 14. This is In reasonable agreement with the results obtained 
In a more rigorous manner In the present study. 

The results presented In Figures 13 and 14 do not differ significantly 
for other r^/r^ ratios. For example, Figures 15 and 16 are contour plots 
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Figure 15. Matrix Stress Contour Plots for rf/rm •• 0.25 
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Figure 15 (continued) 
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Figure 15 (continued) 
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of thr«« dlfforont strooa eoapononts throughout tho alunlnum nuitrlx. 
for rf/rpi ratio* of 0.25 and 0.84, raapaetlvaly, for an appllad axial 
■trait of 4 kal. Thaaa ara computar-ganaratad contour plota; tha 
nunarleal valuta of tha eontoura, and tha Idantlfylng aymbola, llatad In 
tha Itgand banaath aaeh plot, ara uaaful In Idantlfylng closely apacad 
contours. 

Tha octahedral ahaar atraas plots of Figures ISa and 16a ara 
significant in that this strass Is used to daflna yield in the present 
analysis. The contour values have been normalized by dividing by the 
yield value of the octahedral shear stress, viz, 16.97 ksi, as Indicated 
In the caption. As can be seen In both Figures 15a and 16a, the octahedral 
shear stresses, although still In the elastic range (the highest normalized 
values are less than one), are highly localized at the fiber break. It Is 
obvious that first yield will occur at this location. 

The maximum principal stress Is plotted In Figures 15b and 16b. In 
the present examples, the maximum normal (l.e., principal) stress failure 
criterion (see Table 1) was used, as previously noted. The high axial 
normal stress In the aluminum matrix Immediately adjacent to the fiber break 
(previously Indicated In Figure 13 for rf/r^ ■ 0.46) Is obvious, as Is 
the rapid decrease In stress away from this location. 

Shear stress contours In the matrix are shown In Figures 15c and 16c. 
Again, the stress concentration at the fiber break can be readily observed, 
consistent with the plot for rf/r„ - 0.46 presented In Figure 14. 

The stress contour plots of Figures 15 and 16 were for an applied axial 

stress of only 4 ksi, well within the elastic range of the aluminum matrix 

* 

material response (see Figure 4). Figure 17 represents a series of normalized 
octahedral shear stress plots for r£/r„ - 0.46, for increasing levels of 
applied axial stress,” well into the range of Inelastic material response and 




42 


BORON/RLUMINIUM 

NOPNflLIZED OCTR^CORn. SHERR STRESS NORNRLIZED NRT 1697P^I 



INC 1 

RVa STRESS R= 
CONTOUR VRLUES 


TEMPERATURE 0.0 
0.0 4000 

.05 .16 .28 .40 


MOISTURE PERCENT 0.0 


.51 .63 


K M O A * X 

a) Applied Axial Stress Og “ 4,000 psi 


Figure 17. Matrix Octahedral Shear Stress Contour Plots (Normalized) for 
rf/^m “ 0.46, Showing Influence of Crack Propagation 
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Figure 17 (continued). 




BORON/flLUMINIUM 


NORHRLIZEO QCTmCDiVl. SHEAR STRESS NORNflLlZED IRT 1697RSI 


.13 .1313 



INC 6 TEMPERATURE 0.0 MOISTURE PERCENT 0.0 

RVi STRESS Rs 0.0 Z= 12,989 

CONTOUR VRLIES .13 .26 .41 .55 .70 .83 

M K e * X 

c) Applied Axial Stress 5^ •* 12,989 psi 

Figure 17 (continued) 




BORON/FILUMINIUM 


NQRKfLIZED OCTflrtDfiflL SHEAR STRESS NORMflLlEED W!T 1697RSI 


.17 .17.17 



.7Q52 


IMC 8 TErf»ERflTURE 0.0 MOISTU?E PERCENT 

flVG STRESS R= 0.0 Z= 3,696 

CONTOUR VfiLUES .17 .34 .52 .70 .88 1.04 

K K © A * * 

d) Applied Axial Stress = 13,696 psl 


Figure 17 (continued) 
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crack propagation. The normalieed values of the octahedral shear stress 
contours are generally less than one on^y Lecauee these plots represent 
the state of stress immediately following an increment of crack growth 
(during which the stresses decrease, as Indicated in Figure 7). 

Beyond the loading state represented by Figure 17d, the applied 
stress continues to decrease with continued crack propagation (see Figure 7). 

It will be n'^t'cvi in Figure 17 that a high stress concentration persists 
in the region of the propagating crack, as would be expected. 

All of the preceedlng numerical results were generated using the 
maximum normal stress failure criterion. The applied axial stress at 
which the crack initiates, and the ultimate applied stress, are sensitive 
to the failure criterion utilized. Six different failure criteria have been 
Included in the present analysis, as presented In Table 1 of Section 2.3. 
Results using each of these theories, for the boron /aluminum composite 
model with rf/r„, *• 0.46, are presented in Table 5. 

Much more study will be required to understand the full significance 
of the different results obtained. It is obvious, however, that the 
considerable differences exhibited in Table 5 indicate that which failure 
criterion selected is Important, and that it slould be relatively straight- 
forward to select the appropriate one, or ones, by correlation with 
experiment . 

In addition to the variation in predicted stress levels, the 
appropriate failure criterion will also be governed by the correlation 
between predicted and exper ' mentally observed crack propagation patterns. 

For example, the maximum normal stress theory and the octahedral shear 
stress theory, which predict axial tensile ultimate strengths of 21.7 and 
24.2 ksl, respectively (Table 5), both indicate a Mode I (opening) failure 
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T*bl« 5 

Influence of Assumed Feilnre Theory on Crack Propagetion 
In a Boron/Aluminum Composite, r./r_ ■ 0.46 

Z w 


Failure Theory 


.. Maximum Normal Stress 

5,666 

i 

2.1 '96 

*. Maximum Shear Stress 

4,848 

36,038 

1. Octahedral Shear Stress 

7,923 

24,196 

>. Tsal-Hlll 

4,693 

40,356 

>. Hoffman 

2,262 

14,750 

). Tsai-Wu (Modified) 

2,497 

14,225 
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mod*. Th« maxlmus; normal atraaa thaory raaulta waro prai.'antad In Flgura 8. 
Tha oetahadral ahaar atraaa thaory raaulta ara ahown in Flgura 18. Hoda II 
(ahaar) fallura la pradomlnata whan tha maximum ahaar atraaa or Taal-Hill 
failure criterion la aaaumad* tha pradlctad ultimate rtrengtha being 36.0 
and 40.4 kal» reapactlvely. The Hoffman and Taal>Wu failure criteria, 
both of which exhibit predomlnently Model I type fallurea, predicted ultimate 
strengcha of 14.8 and 14.2 kal, reapactlvely. Figure 19 represents the 
crack propagation pattern predicted using the maximum shear stress failure 
criterion. As can be seen, there la a strong tendency for the crack to 
propagate along the fiber-matrix Interface, due to the high shear stress 
In this region (as shown In Figure 14). Of course, for complete failure 
(fracture) of the model composite to occur due to the applied axial 
tensile stress, the matrix crack must eventually propagate to the 
outer radius, as shown In the last sketch of Figure 19. 








21,379 


7,-7,922 7,-24,196 

7,-5,025 7,-17,659 

Figure 18, Crack Propagation In a Boron /Aluminum Composite Model as Predicted 
by the Octahedral Shear Stress Theory, rf/rnt"0.46 (stresses in psl) 



9‘,>4,847 7,*24,200 9,«36,038 

^>4,712 9,«23,232 

Figure 19. Crack Propagation In a Boron/Aluminum Composite Model as Predicted 
by the Maximum Shear Stress Theory, rf/r™ ■ 0.46 (stresses In psl) 
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4.3 Energy of Crack Propagation In a Modal Composite 

An estimation of the absorption of energy associated with crack 
propagation Is an Important aspect of any fracture mechanics study. A 
significant amount of Information can be obtained by studying the variation 
of energy levels during crack propagation 

4.3.1 Evaluation of Fracture Energy 

The total potential «»r. >■ V , of a cracked body consists of elastic- 

strain energy, U , the energy absorbed in the plastic region, U , surface 
e p 

energy, U^, and the potential energy of external loads, A balance of 
these energies during crack propagation, as presented by Griffith, Irwin, 
and Orowan [301 forms the basis of linear elastic fracture mechanics (LEFM), 
The above criterion states that a crack propagates when the energy 
release rate, G, is equal to the energy absorbed, R, i.e., G=R. 

Since the energy release rate and the stress intensity factor, K, 
are related when G=G^ (critical value of G) , the corresponding stress 
intensity factor is K . This critical K value is the so-called fracture 
toughness (i.e,, resistance to cr-^ck growth). A significant amount of 
work has been done to characterize the fracture toughness as a material 
property- It has been observed that the crack resisting parameter R is 
not a constant, but increases as the crack extends due to the Increase in 
the plastic enclave absorbing acre and more energy. The fact that R 
is not a constant, but increa^ae3 as the crack extends, explains the observed 
stable crack growth. Plots of G and R as functions of crack length a (see 
Figure 5), for a range of initial crack lengths, are known as R-curves. 

The R-curves are very useful in understanding the behaviour of crack 




propagation. 
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A. 3.2 Calculation of Total Pofntlal Efttrgy and Energy Reif« Rate 

In the preaent finite element analyala, the total potential energy le 
calculated at each load Increment. The total potential energy, V, Is 
given as 

V - u + n (6) 

where U Is the total strain energy (U plus U_) and (2 Is the potential 
energy of the external loads. Surface energy is assumed negligible. 

The total strain energy U is given by 


U 



Vl 


Uot dVi 


(7) 


where 

^1 


strain energy density In a finite element 
|{o}’’{de) 


n * number of finite elements 
■ volume of an element 


The potential energy of the external loads Is 




( 8 ) 


where 6^ Is the displacement of the loaded boundary point j . 
The energy release rate, G, Is then calculated as 


^1 - ^2 
Aa 


(9) 


where is the total potential energy at the Initiation of the crack 
and V 2 is the total potential energy after the crack has extended an 
Increment Aa. 

Figure 20 shows the variations of total potential energy In the 
boron/aluminum, broken-fiber model (Figure 1) as the crack extends. 


TOTAL POTENTIAL ENER6Y/irrS(in 



Figure 20. Total Potential Energy (Normalized) Versus Crack Extension 
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Slncei in the present finite element analysis » ctacks extend In finite 
lengths (equal to the size of the failed element) » a step variation is 
observed. The dashed lines Indicate the smoothed curves. 

In Figure 21, R-curves corresponding to throe different initial 
crack lengths (equal to the broken fiber radius) are plotted. It can be 
observed that the energy release (i.e., the energy absorbed) is higher 
for the smaller initial crack length, and attains an unstable condition 
earlier than for the larger initial crack length geometries. 


ENERGY RELEASE RATE, G (ENERGY ABSORPTION RATE, R)(IO*in*/in) 




SECTION 5 


DISCUSSION 

The primary purpose in developing the axlsymmetrlc analyals, with 
Inelastic crack propagation capablllty» was to permit the cortelaf .on of 
analytical predictions of crack propagation and failure with exp'^rlmental 
measurements. It Is the experimental work which encourages the use of a 
simple composite model. By using a simple test specimen, experimental 
measurements are relatively easy to perform, and Interpret. 

At the present time, experimental data for the boron/aluminum 
composite, single broken fiber model is not yet available. The original 
plan was that NASA-Lewis would generate this <!ata. It now appears that 
this may not be possible. One alternative Is to undertake this as a 
subsequent study at the University of Wyoming. 

Another alternative Is to model an epoxy-matrix composite, and then 
correlate the predicted response with available single fiber composite 
data. For example, preliminary discussions with Drzal [311 Indicate that 
he may have suitable data for a graphlte/epoxy model composite. To 
establish that the present analysis will perform as well in predicting 
stable crack propagation In an epoxy matrix, a much less ductile material 
than the aluminum matrix used In obtaining the present results, another 
serleF of computer runs were made, for both glass/epoxy and graphlte/epoxy. 
These results are summarized In Reference [32]. In general, the results 
obtained were similar to those presented here. In the sense of stable cracks 
being propagated. One obvious and distinct advantage of using a pol)rmer 
matrix Is Its transparency. The propagating crack can be observed directly, 
and measurements of such parameters as crack opening displacement (the gap 


width between the ends of the broken fiber) can be made optically. For 
the metal matrix composite model, an X-ray technique or something similar 
will be required. 
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